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Microencapsulated  phase  change  materials  (MPCM)  have  been  recognized  as  effective  materials  to 
enhance  heat  transfer,  and  to  improve  heat  storage  performance  in  thermal  energy  system.  Microen¬ 
capsulated  phase  change  material  slurry  (MPCS)  can  apply  to  heat  transport  and  thermal  energy  storage 
systems.  In  order  to  fully  develop  the  application  of  MPCS  in  thermal  energy  system,  more  researches 
on  preparation  and  characteristics  of  MPCS  have  been  done.  This  paper  presents  a  review  on  microen¬ 
capsulation  methods  and  thermal  characteristics  of  MPCS.  It  focuses  on  the  thermal  properties  and  heat 
transfer  characteristics  of  MPCS  flowing  in  horizontal  circular  pipe.  Some  phase  change  materials  and 
microencapsulation  methods  are  analyzed  and  discussed.  Theoretical  models  for  analyzing  heat  trans¬ 
fer  characteristics  of  MPCS  flowing  in  the  pipe  are  presented.  Several  factors  affecting  the  heat  transfer 
characteristics  of  MPCS  are  also  summarized. 
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1.  Introduction 

The  demand  of  energy  becomes  larger  as  the  rapid  development 
of  the  economy  and  human  society.  At  present,  energy  supply  is 
mainly  from  fossil  fuels.  However,  conventional  fossil  fuels  are  lim¬ 
ited,  and  the  use  of  fossil  fuels  results  in  emission  of  greenhouse 
gas  and  many  harmful  gases,  which  cause  the  environmental  pol¬ 
lution  and  climate  changes.  In  recent  decades,  many  researchers 
do  their  efforts  to  deal  with  these  problems.  There  are  three  main 
ways:  (1)  reduce  the  energy  demand,  (2)  improve  the  energy  effi¬ 
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ciency  of  power  systems  and  reduce  the  waste  of  energy,  (3)  find 
new  and  renewable  sources  of  energy.  Thermal  energy  storage  is 
an  effective  method  that  can  reduce  the  mismatch  between  supply 
and  demand,  thus  it  also  improves  the  energy  efficiency  of  power 
systems. 

Thermal  energy  storage  includes  sensible  heat  storage  and 
latent  heat  storage.  Due  to  sensible  heat  storage  material  has  a 
low-energy  storage  density,  it  leads  to  that  heat  storage  system 
requires  large  space  to  store  the  materials,  and  the  construction 
cost  also  increases.  Latent  heat  storage  material  provides  a  high- 
energy  storage  density  when  it  undergoes  a  phase  change  process, 
and  it  is  able  to  store  energy  at  a  constant  temperature  or  with 
a  limited  range  of  temperature  variation,  so  phase  change  mate¬ 
rials  are  widely  used  in  thermal  energy  storage  systems.  But  the 
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Nomenclature 

c  specific  heat  (J /kg  K) 

D  diameter  of  pipe  (m) 

dc  diameter  of  the  core  of  microcapsules  (m) 

dp  diameter  of  microcapsules  (m) 

h  heat  transfer  coefficient  (W/m2  K) 

k  thermal  conductivity  (W/m  K) 

m  mass  flow  rate  (kg/s) 

Nu  Nusselt  number 

Pe  Peclet  number 

Q.  heat  of  fusion  (J/kg) 

q'w  wall  heat  transfer  rate  (W) 

qw  wall  heat  flux  (W/m2 ) 

R  radius  of  pipe 

Re  Reynolds  number 

Rex  local  Reynolds  number 

Ste  Stephan  number 

T  temperature  (°C) 

Th  upper  phase  change  temperature  (°C) 

1/  lower  phase  change  temperature  (°C) 

u  axial  velocity  (m/s) 

v  radial  velocity  (m/s) 

Greek  letters 

a  volume  fraction 

otm  mass  fraction 

rj  viscosity  (Pas) 

p  density  (kg/m3) 

Subscripts 
b  bulk  fluid 

e  effective 

/  carrier  fluid 

i  inlet 

m  mean 

o  outlet 

p  particle 

w  internal  wall  surface 


direct  use  of  phase  change  materials  for  thermal  energy  storage 
application  is  limited  because  the  material  is  prone  to  be  frozen 
on  the  application  wall  when  it  changes  its  phase  from  liquid  to 
solid.  Due  to  the  low  thermal  conductivity  of  most  phase  change 
materials,  the  frozen  layer  can  drastically  decrease  the  efficiency 
of  heat  transfer  [1,2].  In  order  to  solve  the  problem,  suspension  of 
the  phase  change  material  in  a  carrier  liquid  is  developed  by  being 
emulsified  or  microencapsulated  [3]. 

Microencapsulation  is  a  process  by  which  particles  of  solid  or 
liquid  material  are  coated  with  some  polymeric  materials.  The 
capsules  produced  by  microencapsulation  are  in  the  microme¬ 
ter  to  millimeter  range.  The  microencapsulation  technology  has 
been  widely  applied  to  the  pharmaceutical  manufacturing,  chem¬ 
ical  industry  and  biology  engineering.  Microcapsules  with  phase 
change  material  as  core  material  are  mixed  with  some  carrier  fluid, 
so  the  suspension  is  called  microencapsulated  phase  change  mate¬ 
rial  slurry  (MPCS).  The  shell  of  microcapsules  protects  the  core 
material  from  interacting  with  carrier  fluid,  which  enhances  the 
stability  of  the  material  and  avoids  particle  aggregation.  Microen¬ 
capsulation  increases  surface-to-volume  ratio  of  phase  change 
material  and  the  particles  interact  with  particles  or  fluid,  which 
significantly  increases  the  capability  and  heat  transfer  efficiency 
of  MPCS  [3,4].  In  recent  years,  researchers  have  done  much  work 


on  the  preparation  and  properties  of  MPCMS,  and  many  important 
results  were  obtained. 

2.  Preparation  of  microencapsulated  phase  change 
materials 

The  applications  of  microencapsulated  phase  change  material 
(MPCM)  are  mainly  in  thermal  energy  storage,  heat  transfer  and 
temperature  controlling.  For  instance,  it  is  used  in  the  secondary 
loop  of  refrigeration  system  or  applied  to  thermal  insulation  mate¬ 
rial  in  construction.  That  Narita  airport  in  Tokyo  uses  MPCMS  to 
store  cold  energy  is  a  practical  large-scale  application  of  MPCS 
[5].  Griffith  and  Eamesa  [6]  applied  MPCS  in  a  cooled  ceiling  sys¬ 
tem  and  conducted  a  four-month  experiment.  It  was  proved  from 
the  experiment  that  MPCS  could  effectively  reduce  the  flow  rate 
while  keeping  a  constant  temperature.  Mihashi  et  al.  [7]  used 
butyl  stearate  PCM  on  porous  lightweight  aggregates.  The  concrete 
impregnated  could  maintain  interior  temperature  around  18°C. 
Tyagi  et  al.  [8]  and  Zhang  et  al.  [5]  present  reviews  on  the  devel¬ 
opment  of  phase  change  materials  based  on  microencapsulated 
technology.  Microencapsulated  phase  change  materials  are  very 
attractive  in  improving  the  thermal  performance  of  textile  fabrics 
[9,10]. 

Inorganic  phase  change  material  has  some  advantages,  such 
as  large  thermal  conductivity  and  high  energy  storage  density. 
However,  supercooling  is  high  when  the  material  changes  its 
phase  from  liquid  to  solid,  and  phase  separation  happens  in  phase 
change  process.  The  core  material  used  in  MPCM  should  meet  some 
requirements:  (1)  suitable  phase  change  temperature  or  temper¬ 
ature  range,  (2)  large  latent  heat,  (3)  good  thermal  conductivity, 
(4)  low  density  change  and  (5)  low  reactivity.  Paraffin  is  a  kind  of 
practical  organic  material  suitable  for  MPCM.  Mixture  of  organic  or 
inorganic  PCM  can  be  used  to  develop  ideal  core  material. 

There  are  many  kinds  of  organic  or  inorganic  materials  which 
can  be  used  for  the  shell  material.  However,  polymers  are  com¬ 
monly  used.  The  shell  material  should  be  selected  according  to 
the  physical  properties  of  the  core  material.  If  the  core  mate¬ 
rial  is  hydrophilic,  hydrophobic  polymer  is  generally  used  as  the 
shell  material,  and  vice  versa.  The  shell  material  should  meet 
some  requirements  [1 1  ]:  (1 )  good  flexibility,  (2)  good  sealing  tight¬ 
ness,  (3)  endurance  and  (4)  low  reactivity.  Some  polymers  are 
often  selected  as  shell  material,  such  as  polystyrene,  polymethyl¬ 
methacrylate,  Arabic  gum,  gelatin,  amino  plastics,  gelatin-gum 
Arabic,  urea  formaldehyde  resin,  melamine  formaldehyde  resin, 
gelatin  formaldehyde  resin,  and  so  on  [9,1 2-1 6].  The  silica  prepared 
by  hydrolysis  condensation  reaction  of  TEOS  is  also  used  as  shell 
material  [17]. 

Microcapsules  can  be  achieved  by  many  methods,  which  include 
physical  method,  chemical  method  and  physicochemical  method. 
Microencapsulation  methods  of  phase  change  materials  are  pre¬ 
sented  in  Table  1.  The  diameter  ranges  of  microcapsules  are  shown 
in  Table  2  [31]. 

•  Physical  methods :  Spray  drying,  spray  cooling,  air  suspension 
coating,  supercritical  fluid  method,  centrifugal  extrusion,  elec¬ 
trostatic  precipitation,  and  so  on. 

•  Chemical  methods :  Interfacial  polymerization,  in  situ  polymeriza¬ 
tion,  piercing-solidifying,  suspension  cross-link  method,  and  so 
on. 

•  Physicochemical  method :  Phase  separation. 

Fig.  1  presents  type  of  microcapsules,  such  as  mononuclear 
microcapsules:  a  single  core  wrapped  with  a  continuous  shell 
material,  poly-nuclear  microcapsules:  many  cores  coated  with 
a  continuous  shell  material  and  multi-film  microcapsules:  a 
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Table  1 

Microencapsulation  methods  of  phase  change  materials. 


References 

Core  material 

Shell  material 

Microencapsulation  method 

Zou  et  al.  [18] 

Hexadecane 

TDI  and  EDA 

Interfacial  polymerization 

Fang  et  al.  [17] 

Paraffin  wax 

Si02 

Sol-gel  method 

Chen  et  al.  [19] 

1  -Bromohexadecane 

Amino  plastics 

Ozonur  et  al.  [20] 

Coco-fatty  acid 

Gelatin-gum  Arabic 

Coacervation 

Alvarado  et  al.  [21  ] 

Tetradecane 

Gelatin 

Coacervation 

Choi  et  al.  [22] 

Tetradecane 

Melamine  formaldehyde 

In  situ  polymerization 

Alkan  et  al.  [23] 

Docasane 

Polymethylmethacrylate 

Emulsion  polymerization 

Alkan  et  al.  [24] 

Eicosane 

Polymethylmethacrylate 

Emulsion  polymerization 

Hawlader  et  al.  [25] 

Paraffin  wax 

Gelatin  acacia 

Coacervation/spraying  drying 

Fang  et  al.  [26] 

Tetradecane 

Urea  formaldehyde 

In  situ  polymerization 

Zhang  et  al.  [27] 

Octadecane 

Melamine  formaldehyde 

In  situ  polymerization 

Bayes-Garcia  et  al.  [28] 

Mixture  of  alkanes 

Gelatin-gum  Arabic/sterilized  gelatin-gum  Arabic 

Coacervation 

Ai  et  al.  [29] 

Hexadecane 

Polystyrene 

Suspension  polymerization 

Onder  et  al.  [30] 

Hexadecane 

Gelatin-gum  Arabic 

Complex  coacervation 

Octadecane 

Nonadecane 

Fig.  1.  Type  of  microcapsules,  (a)  Mononuclear  microcapsule,  (b)  poly-nuclear  microcapsule  and  (c)  multi-film  microcapsule. 


Fig.  2.  Microscope  profiles  of  MPCM  made  by  different  methods,  (a)  Spray  drying  method  [25],  (b)  coacervation  method  [25]  and  (c)  in  situ  polymerization  method  [26]. 
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Table  2 

The  diameter  range  of  microcapsules  [31  ]. 


Method  of  microencapsulation 

Core  material 

Diameter  range  (p,m) 

Phase  separation 

Solid,  liquid 

2-1200 

Interfacial  polymerization 

Solid,  liquid 

2-2000 

Spray  drying/cooling 

Solid,  liquid 

6-600 

Centrifugal  extrusion 

Solid,  liquid 

1-5000 

In  situ  polymerization 

Solid,  liquid 

1-2000 

Suspend  wrapping  in  air 

Solid 

35-5000 

Electrostatic  precipitation 

Solid,  liquid 

1-50 

Suspension  cross-link  method 

Solid,  liquid 

2-4000 

continuous  core  coated  with  multilayer  continuous  shell  material. 
The  microscope  profiles  of  MPCM  made  by  different  methods  are 
shown  in  Fig.  2. 


3.  Thermal  properties  of  MPCS 


The  thermal  properties  of  MPCS  are  different  from  the  PCM 
and  the  carrier  fluid.  In  general,  several  characteristics  of  MPCS 
are  important  for  design  of  MPCS  system,  such  as  density,  thermal 
conductivity,  specific  heat  capacity  and  viscosity. 

According  to  conservation  of  mass: 


1 

Otm,c / Pc  +  (1  —  Otm,c) / Ps 

(i) 

1 

Otm,p/ Pp  +  (1  —  Ctm,p)/ Pf 

(2) 

The  densities  of  microcapsule  and  carrier  fluid  are  similar  to 
ensure  the  long-term  stability  of  MPCS  while  preparing  the  sus¬ 
pension,  so  am^a.  Though  the  density  of  core  material  varies  by 
10-15%  as  phase  change  process  happens,  the  density  of  slurry  can 
be  analyzed  as  a  constant  because  the  variation  is  less  than  1-2%  in 
low  volumetric  concentration  [32]. 

The  specific  heat  capacity  should  be  evaluated  in  order  to  ana¬ 
lyze  the  phase  change  effects. 

According  to  conservation  of  energy  in  solid  phase  and  liquid 


phase  of  PCM: 

Cp  =  Otm,cCc  +  (1  —  Otm^c)Cs 

(3) 

Cb  =  C%m,pCp  +  (1  —  Otm^p)Cf 

(4) 

Though  the  specific  heat  capacity  is  different  between  solid 
phase  and  liquid  phase,  it  can  be  analyzed  as  a  constant  because 
of  the  low  concentration  of  PCM. 

In  the  theoretical  model,  the  phase  change  material  is  assumed 
to  melt  over  a  temperature  range.  The  specific  heat  capacity  in  the 
melting  process  is  defined  as  the  following  equation: 

fh 

Q=  /  cdT  (5) 

JTi 

Some  works  show  that  the  shape  of  specific  heat 
capacity-temperature  curve  has  a  little  effect  on  the  heat 
transfer  process  [33,34].  Hu  and  Zhang  [35]  analyzed  four  different 
capacity-temperature  curve  (left  triangle,  right  triangle,  rectan¬ 
gular  and  sine  curve),  and  the  result  indicates  that  four  different 
specific  heat  capacity-temperature  curves  have  little  differences 
in  heat  transfer  investigation  except  in  the  thermal  entry  region. 
So  the  specific  heat  capacity  inside  the  melting  temperature  range 
can  be  given  by  the  following  equation: 


Ce  —  Cjj  + 


d 


Therefore  the  specific  heat  capacity  in  the  whole  process  of  heat 
transfer  is  shown  as  follows: 

/ 


1 


Ce  —  Cjj 
Ce  =  Cjj  + 
Ce  —  Cfj 


Q 


h 


Tj  <  T  <  T, 
Ti<T<Th 
Th  <  T  <  T0 


(7) 


The  thermal  conductivity  of  the  microcapsules  is  calculated  by 
the  composite  sphere  approach  [36]: 


1 


kpdp 


1  dp -d 
+  p 


kc  dc 


ksdpdc 


dp  \  —  P c ^  — 


dc 


PsCtm,c 


(8) 


(9) 


The  thermal  conductivity  of  the  MPCS  is  calculated  by  Maxwell’s 
relation  [37]: 


2  +  kp/kf  +  2a(kp/kf  —  1) 
<b  =  f  2  +  kp/k/-a(kp/k/-l) 


(10) 


However,  the  thermal  conductivity  predicted  by  Maxwell’s  rela¬ 
tion  is  lower  than  the  effective  thermal  conductivity  on  flowing 
condition  because  of  the  interaction  between  the  particles  and  fluid. 
The  enhancement  of  the  effective  thermal  conductivity  is  relative 
to  the  particle  size,  shear  rate,  thermal  diffusivity,  concentration 
of  the  slurry,  and  so  on.  The  effective  thermal  conductivity  can  be 
written  as  follows  [38]: 


ke 


=  h»S 


f  =  1  +  BcePe™  =  1  +  Ba  8m 


(11) 

(12) 


{£  =  3.0  m  =  1 .5  Pep  <  0.67 

£  =  1.8  m  =  0.18  0.67  <  Pep  <  250  ^3^ 

£  =  3.0  m  =  ^  Pep  >  250 

The  viscosity  of  MPCS  is  important  to  analyze  the  flowing  pro¬ 
cess  and  the  pressure  drop.  The  pumping  power  increases  as  the 
pressure  drop  increases,  which  is  detrimental  to  the  practical  appli¬ 
cation  of  MPCS.  Higher  viscosity  means  higher  pumping  power  and 
lower  turbulence  of  the  slurry,  so  the  heat  transfer  coefficient  may 
decrease  beyond  the  phase  change  process.  The  viscosity  of  MPCS 
is  affected  by  several  factors,  such  as  the  viscosity  of  carrier  fluid, 
the  microcapsules  concentration,  the  particle  size  and  the  surface 
roughness  of  the  microcapsules.  Some  results  [39,40]  indicated  that 
the  viscosity  was  1 .2-1 1  times  higher  than  that  of  water  while  the 
particle  concentration  increased  from  5%  to  30%.  The  MPCS  can  usu¬ 
ally  be  treated  as  homogeneous  fluid  when  the  volume  fraction  is 
not  high  (up  to  37%).  The  viscosity  can  be  given  as  follows  [41  ]: 


T)b  2.5 

—  =  ( 1  —  a  -  Acr ) 
df 


(14) 


where  A  is  a  parameter  which  is  relative  to  the  size,  the  shape, 
the  rigidity  and  the  type  of  the  microcapsules.  Mulligan  et  al.  [42] 
got  the  value  A  =  3.4  for  the  MPCS  with  particle  of  10-30  p^m  in 
diameter.  Wang  et  al.  [40]  estimated  A  =  4.45  while  investigating 
the  microencapsulated  1-bromohexadecane  slurry  with  particle  of 
10.112  |jim  in  average  diameter  at  different  concentrations.  Yam- 
agishi  et  al.  [39]  analyzed  the  microencapsulated  octadecane  slurry 
with  particle  of  6.3  p>m  in  average  diameter,  and  obtained  A  =  3.7. 
Charunyaorn  et  al.  [38]  indicated  the  value  A  =  1.1 6  when  the  vol¬ 
ume  concentration  was  lower  than  20%. 

The  MPCS  can  be  treated  as  Newtonian  fluid  when  the  volume 
fraction  is  lower  than  30%  [39].  Yamagishi  et  al.  [43]  carried  out 
tests  of  the  MPCS  using  the  microcapsules  with  the  core  mate¬ 
rial  of  n-tetradecane  and  n-dodecane.  As  shown  in  Fig.  3  [39],  the 


Th  -  T, 
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Reb  [.] 

Fig.  3.  The  relation  between  friction  factor  and  the  Reynolds  number  of  MPCS  [39]. 


relation  between  friction  factor  and  the  Reynolds  number  of  MPCS 
is  presented,  and  the  slurry  evidently  displays  the  similar  char¬ 
acteristics  of  Newtonian  fluid.  Yamagishi  et  al.  [44]  indicated  that 
anionic  surfactants  could  decrease  the  viscosity  of  MPCS  effectively 
at  high  mass  fraction  and  suggested  that  the  shape  and  rigidity 
of  the  microcapsules  may  have  an  effect  on  the  viscosity.  Wang 
et  al.  [40]  developed  experiment  to  study  the  rheological  behav¬ 
iors  of  MPCS,  and  the  result  showed  that  the  shear  stress  increased 
with  shear  rate  linearly  with  mass  fraction  below  27.6%.  It  is  also 
found  that  the  states  of  phase  change  material  have  nearly  no  effect 
on  the  rheological  behaviors,  because  the  particles  contacting  with 
the  carrier  fluid  were  always  solid  state.  So  the  viscosity  is  almost 
independent  of  temperature  [21,40,43].  Zhang  et  al.  [27]  carried 
out  comprehensive  studies  on  the  fabrication  and  properties  of 
microcapsules  containing  n-octadecane.  The  result  showed  that  the 
particle  size  and  morphology  were  significantly  relative  to  the  stir¬ 
ring  rate.  As  the  stirring  rate  increases,  the  average  diameter  of 
particles  decreases  and  the  surface  of  the  particle  becomes  smooth. 

For  the  long-term  application,  the  MPCS  should  be  endurable. 
Ohtsubo  et  al.  [45]  operated  several  experiments  to  study  why  the 
microcapsules  break,  and  the  results  presented  that  the  breakage 
of  microcapsules  increased  while  the  ratio  of  particle  diameter  to 
thickness  increased.  Yamagishi  et  al.  [43  ]  indicated  that  the  smaller 
particle  was  more  endurable  by  testing  the  breakage  of  particles 
with  four  different  size.  Alvarado  et  al.  [21]  studied  the  microcap¬ 
sules  made  by  microencapsulating  n-tetradecane  with  gelatin,  and 
indicated  that  the  size  should  be  less  than  1 0  |jim  to  ensure  the  dura¬ 
bility  of  the  particles.  The  experiment  also  showed  the  durability 
increased  as  particle  size  decreased.  Roy  and  Sengupta  [46]  ana¬ 
lyzed  the  stability  of  some  microcapsules  of  eicosane  in  the  range 
of  100-250  |jim,  and  the  result  indicated  that  the  smaller,  thicker- 
walled  particles  were  more  stable  than  the  larger,  thinner  ones.  Ai 
et  al.  [29]  developed  a  method  using  casein  to  protect  the  microcap¬ 
sules  from  breaking  in  the  fabrication  process,  but  the  result  also 
showed  that  the  melting  temperature  increased  and  the  latent  heat 
decreased.  Many  studies  have  focused  on  the  thermal  and  chemical 
stability  of  the  MPCS  when  the  MPCS  is  applied  practically.  Alkan 
et  al.  [23]  carried  out  experiments  to  analyze  the  stability  of  MPCS 
with  n-docasane  as  core  material.  The  result  showed  that  there  was 
no  significantly  change  in  the  properties  of  MPCS  after  5000  ther¬ 
mal  cycles.  Ozonur  et  al.  [20]  used  coco  fatty  acid  mixture  as  core 
material  to  prepare  MPCS.  It  is  proved  that  the  chemical  structure 
of  the  phase  change  material  does  not  change  and  the  geometric 
profile  of  the  core  material  remains  the  same  after  50  cycles.  Some 
other  papers  have  also  been  published  to  analyze  the  thermal  and 
chemical  stability  of  MPCS  [18,25,47,48]. 


4.  Heat  transfer  model  of  MPCS  in  the  circular  pipe 


The  specific  heat  capacity  is  assumed  to  be  a  function  of  the 
temperature,  so  the  governing  equation  can  be  given  as  follows: 


dT  dT 

U  —  +V^~ 

dz  dr 


1  d 

rpbcb  dr 


(15) 


The  energy  equations  that  govern  heat  transfer  in  system  are 
based  on  the  following  assumptions:  (1)  the  slurry  is  a  uniform 
Newtonian  fluid,  and  the  physical  properties  are  constant  except 
the  specific  heat  capacity,  (2)  the  flow  is  fully  developed  laminar 
flow,  so  the  radial  flow  can  be  neglected,  (3)  the  interfacial  ther¬ 
mal  resistance  of  the  microcapsules  can  be  neglected,  (4)  the  inlet 
temperature  of  the  slurry  is  lower  than  the  upper  phase  change 
temperature,  (5)  the  axial  conduction  and  viscous  dissipation  can 
be  neglected.  So,  the  formulation  can  be  written  as: 

dr  i  d 

U  dz~  rpbcb  dr 


The  boundary  conditions  are  given  as  follows: 


lz=0 


=  Ti 


dT 

dr 


=  0 


r= 0 


dT 

dr 


Qw 


r=R 
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(17) 


e.w 


Define  the  dimensionless  temperature,  radius  and  length  of  the 
pipe,  respectively: 


e  =  T-Tj  f  =  r 
qwR/kb  R 


z  = 


R 


(18) 


The  velocity  profile  of  the  laminar  flow  is  given  as: 
u  =  2um(l  -  r'2) 
m 

Urn  — 


mlpb 


(19) 

(20) 


The  dimensionless  governing  equation  can  be  written  as: 

do 


(  ,2^30  19 

Peh(  1  -  r  )  —  =  — — 
b  Jdz '  r'  dr' 


r'f 


dr' 


(21) 


Peb  =  RebPrb  = 


Dumpb  cb 
kb 


(22) 


The  dimensionless  boundary  conditions  are  given  as  follows: 


do_ 

dr' 


=  0 

r'= 0 


do 


(23) 


In  the  heat  transfer  process  of  MPCS  in  the  circular  pipe, 
the  phase  change  material  goes  through  three  phases:  solid, 
solid-liquid  and  liquid,  respectively.  So  the  local  mean  temperature 
of  the  MPCS  in  the  whole  process  is  not  linear.  Choi  et  al.  [49]  pro¬ 
posed  a  “three-region  melting  model”  to  calculate  the  local  mean 
temperature. 


Region  1  (solid  region)  :  L\  =  ™Cb^1 — Inhlli 

Qw 

Region  2  (phase  change  region)  :  L2  =  L-L\  -  L3 

Region 3  (liquid  region)  :  L\  =  — Ihli 

qw 


The  local  mean  temperature  can  be  given  as: 

z  0  <  z  <  L\ 


'  j  _T  ,^1  Tmi 

1  m  —  1  i  i-  ■ 


< 


Tm  —  Tj  + 


u 

Th-T, 
^2 


(z-L,) 


Lj  <  z  <  Lj  +  L2 


Tm  —  Tb  + 


Tm,  o  —  Tf j 
i-3 


(z-Li-L2)  L\  +L2  <z  <L 


(24a) 

(24b) 

(24c) 


(25) 
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Fig.  4.  The  variation  of  dimensionless  wall  temperature  with  dimensionless  axial 
length  [19]. 


To  measure  the  heat  transfer  efficiency  of  MPCS,  several  param¬ 
eters  are  used: 


Heat  transfer  coefficient :  hz 
Local  Nusselt  number  :  Nuz 


hzD 

kb 


(26) 

(27) 


Dimensionless  local  temperature  :  0Z 


Tz  —  Tj  q  —  ~  ^i 

qwR/kb  w  ~  qwR/kb 

(28) 


5.  Effect  of  factors  on  heat  transfer  characteristics  of  MPCS 

MPCS  is  advantageous  compared  to  the  conventional  sensible 
heat  transfer  fluid  because  of  the  latent  heat  effect  of  the  phase 
change  material.  The  specific  heat  capacity  is  much  higher,  which 
can  greatly  improve  the  heat  transfer  efficiency  of  the  slurry. 

In  the  previous  researches,  some  numerical  and  experimental 
analyses  were  made  to  estimate  the  performances  of  MPCS.  Kasza 
and  Chen  [3]  indicated  that  the  heat  transfer  efficiency  of  MPCS 
could  be  improved  to  3  times  than  that  of  the  conventional  fluid  by 
studying  the  performance  of  solar  energy  or  waste  heat  utilization 
system  with  phase  change  slurry.  Charunyaorn  et  al.  [38]  carried 
out  a  numerical  study  of  the  MPCS  flowing  in  the  pipe,  and  the 
result  estimated  that  the  improvement  of  the  heat  transfer  coeffi¬ 
cient  can  be  4  times  higher  than  that  of  single  phase  fluid.  Goel  et  al. 
[50]  conducted  experiment  to  analyze  the  laminar  forced  convec¬ 
tion  heat  transfer  of  MPCS,  and  the  results  were  compared  with  that 
of  Charunyaorn  et  al.  [38].  Though  the  results  only  agreed  qualita¬ 
tively  with  the  numerical  ones  and  the  difference  was  about  45%, 
those  indicated  that  the  dimensionless  wall  temperature  could  be 
50%  reduction  to  that  of  single  phase  fluid.  Zhang  and  Faghri  [51] 
analyzed  the  results  of  Charunyaorn  et  al.  [38]  and  Goel  et  al.  [50] 
and  pointed  out  the  supercooling  was  the  most  important  reason  to 
result  in  the  difference.  Chen  et  al.  [  1 9]  developed  an  experiment  to 
analyze  the  performance  of  the  MPCS  formed  by  microencapsulat¬ 
ing  industrial  grade  1-bromohexadecane  as  phase  change  material. 
As  shown  in  Fig.  4,  the  decrease  of  the  dimensionless  internal  wall 
temperature  for  15.8wt%  MPCS  could  be  30%  higher  than  that  of 
water.  Choi  et  al.  [49]  studied  the  phase  change  emulsion  in  turbu¬ 
lent  flow  under  constant  heat  flux.  The  result  indicated  that  heat 
transfer  coefficient  of  the  emulsion  increased  in  region  1 ,  decreased 
in  region  2  and  increased  again  in  region  3.  Several  experiments  of 
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Fig.  5.  Heat  transfer  coefficient  with  distance  [21  ]. 


thermal  performance  of  MPCS  showed  similar  results  [21,39,40], 
which  is  presented  in  Fig.  5  [21  ].  Alvarado  et  al.  [21  ]  estimated  that 
the  heat  transfer  coefficient  reached  the  maximum  near  the  melting 
point.  Wang  et  al.  [52]  carried  out  experiments  to  investigate  the 
heat  transfer  performances  of  MPCS  in  a  circular  horizontal  tube 
with  constant  heat  flux  at  both  laminar  and  turbulent  flow  con¬ 
ditions.  It  is  found  that  the  heat  transfer  coefficients  for  laminar 
flow  are  always  higher  than  that  of  single  phase  fluid.  The  average 
Nusselt  number  of  turbulent  slurry  can  be  2.5  times  higher  than 
that  of  water.  Two  heat  transfer  correlations  were  also  proposed 
based  on  the  experimental  data  to  describe  the  performances  of 
MPCS  at  laminar  flow  and  turbulent  flow  conditions,  respectively. 
Ozonur  et  al.  [20]  designed  experiments  to  study  the  heat  transfer 
rate  of  the  microencapsulated  phase  change  material  and  the  pure 
phase  change  material  (coco  fatty  acid).  The  result  showed  that 
the  heat  transfer  rate  of  microencapsulated  phase  change  material 
was  higher  compared  to  that  of  pure  phase  change  material  during 
charging  process. 

Hu  and  Zhang  [35]  presented  a  numerical  analysis  of  convec¬ 
tive  heat  transfer  enhancement  with  MPCS.  It  is  modified  to  define 
the  local  Nusselt  number  because  the  conventional  local  Nusselt 
number  which  is  used  to  describe  the  single  phase  fluid  cannot 
accurately  describe  the  performance  of  MPCS.  The  energy  equation 
is  analyzed  and  the  result  shows  that  the  modified  local  Nusselt 
number  is  a  function  of  several  parameters,  such  as  Re,  Pr,  spe¬ 
cific  heat,  the  uniformity  of  velocity  profile  and  the  dimensionless 
temperature  gradient.  But  the  parameters  are  not  independent  and 
some  are  strongly  relative  to  the  properties  of  microcapsules  and 
MPCS.  However,  they  made  numerical  simulations  and  predicted 
that  Ste  number  and  mass  fraction  were  the  most  important  param¬ 
eters  affecting  the  heat  transfer  enhancement.  The  supercooling, 
phase  change  temperature  range,  the  particle  size  and  the  Re  num¬ 
ber  have  influence  on  the  heat  transfer  performance.  Zhao  et  al. 
[53]  made  parametric  analysis  of  enhanced  heat  transfer  for  lami¬ 
nar  flow  of  MPCS  in  a  circular  tube  with  constant  wall  temperature. 
It  indicated  that  the  result  is  similar  to  that  of  Hu  and  Zhang  [35]. 
Six  parameters  such  as  Ste  number,  volume  fraction,  supercool¬ 
ing,  phase  change  temperature  range,  ration  of  particle  size  to  tube 
radius  and  Re  number  are  relative  to  the  performance  of  MPCS, 
and  the  Ste  number,  volume  fraction  and  particle  size  are  the  most 
important. 
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Fig.  6.  Effect  of  Ste  on  dimensionless  wall  temperature  [50]. 


5.2.  Effect  of  Ste  number  on  heat  transfer  characteristics  ofMPCS 


The  Ste  number  is  a  dimensionless  number  to  present  the  ratio  of 
the  sensible  heat  to  the  latent  heat,  which  is  given  as  Ste  =  ce{qwr/kb ). 
The  influence  caused  by  Ste  is  strongly  dependent  of  the  phase 
change  process.  Charunyaorn  et  al.  [38]  carried  out  a  numerical 
study  of  the  MPCS  and  pointed  out  that  the  Ste  number  was  the  most 
important  parameter  affecting  the  heat  transfer  process.  Goel  et  al. 
[50]  carried  out  experiment  under  laminar  flow  and  constant  heat 
flux,  and  the  data  also  showed  that  Ste  number  was  the  most  impor¬ 
tant  parameter,  especially  for  Ste  less  than  1.0,  as  shown  in  Fig.  6. 
Zeng  et  al.  [54]  developed  an  experiment  to  study  the  heat  trans¬ 
fer  characteristics  ofMPCS  in  laminar  flow,  and  the  result  showed 
Ste  number  was  the  most  dominant  parameter  to  affect  the  local 
Nusselt  number.  Zhang  et  al.  [55]  carried  out  theoretical  analysis  of 
convective  heat  transfer  enhancement  ofMPCS.  The  numerical  sim¬ 
ulation  results  presented  that  the  fluctuation  amplitude  and  range 
increase  as  the  Ste  number  decreases  in  the  region  2.  Chen  et  al. 
[19]  indicated  that  the  expression  Ste  =  ce(qwr/kb)  could  not  show 
the  sensible  heat  in  the  heat  transfer  process,  so  a  new  Ste  number 
was  defined  as: 


Ste  = 


Ce(T0  —  Tj)  —  Q 

a 


Qi 


w 


liiQ 


-  1 


(29) 


5.2.  Effect  of  mass  fraction  on  heat  transfer  characteristics  of 
MPCS 


The  mass  fraction  affects  the  thermal  conductivity  and  the  vis¬ 
cosity  ofMPCS  according  to  Eq.  (10).  Mulligan  et  al.  [42]  carried  out 
experiments  to  study  the  properties  of  MPCS  with  n-octadecane, 
n-eicosane,  n-heptadecane  and  n-dodecane  as  core  material.  It  is 
found  the  diameter  of  microcapsules  is  between  10  and  30  p>m  and 
the  specific  heat  increases  significantly  as  mass  fraction  increases. 
Yamagishi  et  al.  [43]  conducted  several  experiments  under  turbu¬ 
lent  flow  with  MPCS  of  which  the  size  was  from  2  to  10  p>m.  The 
results  showed  that  the  performance  of  the  slurry  changed  from 
turbulent  flow  to  laminar  flow  at  high  mass  fraction  and  the  pres¬ 
sure  drop  decreased.  Experimental  data  of  Yamagishi  et  al.  [39] 
shows  that  heat  transfer  coefficient  increases  with  mass  fraction 
at  the  same  Re  number  (as  shown  in  Fig.  7),  and  decreases  with 
higher  heating  rate.  The  experimental  results  of  Inaba  et  al.  [57] 
indicated  that  the  ratio  between  heat  transported  and  pumping 
power  decreased  while  increasing  the  mass  fraction  of  large  sized 
particles  at  laminar  flow  condition. 

5.3.  Effect  of  Re  number  on  heat  transfer  characteristics  ofMPCS 

Re  number  presents  the  turbulence  ofMPCS  which  can  improve 
or  suppress  the  heat  transfer  coefficient.  Roy  and  Avanic  [56]  exper¬ 
imental  results  indicate  the  Re  number  has  a  significant  effect  on 


test  section  (X  /  D)  [-] 

Fig.  7.  Local  heat  transfer  coefficients  of  MCPCM  slurries  for  different  Reynolds 
numbers  [39].  (a)  Rez  =  931 1-1 1,020;  c  =  0.12,  (b)  Rez  =  11,094-12,962;  c  =  0.07,  (c) 
Rez  =  5294-6329;  c  =  0.1 5  and  (d)  Rez  =  5488-6495;  c  =  0.12. 


heat  transfer  performance  with  phase  change  material  emulsions. 
Inaba  et  al.  [57]  conducted  experiments  to  study  the  MPCS  with 
large  sized  particles  of  high  mass  concentration.  The  results  showed 
that  the  pressure  drop  decreased  at  turbulent  flow  and  increased 
at  laminar  flow.  Numerical  analysis  of  convective  heat  transfer 
enhancement  with  MPCS  at  laminar  flow  condition  of  Hu  and  Zhang 
[35]  indicated  that  the  average  Nu  number  increased  with  increas¬ 
ing  Re  for  the  same  heat  transfer  fluid.  Zeng  et  al.  [54]  indicates 
that  the  dimensionless  wall  temperature  decreases  as  Re  increases 
at  the  same  dimensionless  length. 


5.4.  Effect  of  supercooling  and  phase  change  temperature  range 
on  heat  transfer  characteristics  ofMPCS 

Supercooling  is  the  temperature  difference  between  freezing 
and  melting  points.  Researchers  explain  the  supercooling  with  the 
classical  nucleation  theory  which  indicates  that  the  phase  change 
process  happens  because  of  the  homogeneous  or  heterogeneous 
nucleation  mechanism.  Montenegro  et  al.  [58]  found  that  homoge¬ 
neous  nucleation  mechanism  caused  greater  supercooling.  Some 
experimental  results  show  that  heterogeneous  nucleation  is  pre¬ 
ferred  to  avoid  undesirable  supercooling  [43,59].  The  numerical 
analysis  of  Hu  and  Zhang  [35]  estimates  that  the  heat  transfer 
enhancement  increases  as  the  supercooling  and  phase  change 
temperature  range  decrease.  Yamagishi  et  al.  [43]  carried  out 
experiments  and  indicated  that  supercooling  and  the  phase  change 
temperature  range  increased  as  the  particle  size  was  less  than 
100  |xm,  and  it  was  presented  that  using  nucleating  agents  with 
molecular  structure  similar  to  the  core  material  could  considerably 
suppress  the  supercooling.  Alvarado  et  al.  [21  ]  prepared  MPCS  with 
microcapsules  containing  94%  tetradecane  and  6%  tetradecanol, 
and  the  experimental  data  indicated  that  the  supercooling  signif¬ 
icantly  decreased.  However,  the  latent  heat  of  fusion  decreases 
significantly  as  nucleating  agent  concentration  increases  [60].  Fan 
et  al.  [61]  conducted  experiment  to  study  the  properties  ofMPCS 
with  octadecane  as  core  material  and  sodium  chloride,  octadecanol 
and  paraffin  as  nucleating  agents,  respectively.  The  results  show 
that  the  supercooling  decreases  significantly  when  the  nucleating 
agent  increases  at  proper  concentration.  Zhang  et  al.  [55]  indicates 
that  the  fluctuation  amplitude  increases  with  decreasing  phase 
change  temperature  range  when  phase  change  happens. 
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Fig.  8.  Dimensionless  wall  temperature  for  different  particle  size  [54]. 


5.5.  Effect  of  particle  size  on  heat  transfer  characteristics  ofMPCS 

The  particle  size  is  also  relative  to  the  thermal  conductivity 
of  MPCS  according  to  Eq.  (8).  In  previous  researches,  Hetsroni 

[62]  indicated  that  the  particles  size  affected  the  heat  transfer  by 
enhancing  or  suppressing  the  turbulence  of  the  slurry.  Liu  et  al. 

[63]  analyzed  the  relation  of  heat  transfer  coefficient  with  the  size 
of  particles  added  into  single  phase  fluid.  The  results  estimated  that 
there  was  a  critical  size,  and  when  the  particles  were  larger  than  the 
critical  size,  the  heat  transfer  coefficient  would  be  enhanced,  and 
vice  versa.  Roy  and  Avanic  [56]  carried  out  an  experiment  to  study 
the  laminar  heat  transfer  with  phase  change  material  emulsion. 
The  results  show  the  characteristics  are  similar  to  MPCS,  suggest¬ 
ing  that  the  shell  of  particle  seems  to  have  no  significant  effect  on 
heat  transfer  process.  Inaba  et  al.  [57]  analyzed  MPCS  with  different 
sized  particles.  The  results  indicate  that  the  heat  transfer  perfor¬ 
mance  ofMPCS  with  different  size  is  better  than  that  of  single  sized 
ones.  The  experimental  results  of  Zhang  et  al.  [27]  presented  that 
the  particle  size  had  no  effect  on  the  melting  behaviors  of  microcap¬ 
sules,  but  significantly  affected  the  freezing  behaviors.  Zeng  et  al. 
[54]  results  showed  that  the  effect  caused  by  particle  size  was  inde¬ 
pendent  of  phase  change,  and  the  dimensionless  wall  temperature 
decreased  as  the  particle  size  increased,  which  is  shown  in  Fig.  8. 

6.  Conclusions 

A  review  of  MPCS  used  for  heat  transfer  in  a  circular  tube  has 
been  carried  out  in  this  paper.  The  materials,  fabrication  and  prop¬ 
erties  of  MPCS  are  presented.  The  microencapsulation  technology 
is  very  useful  to  apply  phase  change  materials  into  thermal  stor¬ 
age  and  heat  transfer  field.  The  heat  transfer  characteristics  and 
parameters  affecting  the  heat  transfer  process  are  shown,  which 
provides  some  basic  information  on  the  research  of  MPCS.  The 
previous  researches  prove  that  the  MPCS  makes  the  heat  trans¬ 
fer  coefficient  improving.  However,  more  work  need  to  be  done  for 
developing  the  practical  application  ofMPCS. 
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